Draft version November 6, 2012 

Preprint typeset using L^T^X style emulatcapj v. 5/2/11 



(N 

o 

(N 

> 
O 



o 

u 

6 



> 
00 

o 



X 



SUPERMASSIVE SEEDS FOR SUPERMASSIVE BLACK HOLES 
Jarrett L. Johnson 1 , Daniel J. Whalen 2 , Hui Li 1 and Daniel E. Holz ,! 

1 Nuclear and Particle Physics, Astrophysics and Cosmology Group (T-2), 
Los Alamos National Laboratory, Los Alamos, NM 87545; jlj@lanl.gov 
2 McWilliams Fellow, Department of Physics, Carnegie Mellon University, Pittsburgh, PA 15213 and 
3 Enrico Fermi Institute, Department of Physics, and Kavli Institute for Cosmological Physics, 
University of Chicago, Chicago, IL 60637 
Draft version November 6, 2012 

ABSTRACT 

Recent observations of quasars powered by supermassive black holes (SMBHs) out to z > 7 allow 
to constrain both the initial seed masses and the growth of the most massive black holes (BHs) in the 
early universe. The combination of the limited role of mergers in growing seed BHs as inferred from 
recent cosmological simulations, the sub-Eddington accretion rates of BHs expected at the earliest 
times, and the large radiative efficiencies of the most massive BHs inferred from observations of 
active galactic nuclei at high redshift, all suggest that the initial BH seeds may have been as massive 
as > lO 5 Af . This is consistent with the prediction of the direct collapse scenario of SMBH seed 
formation, in which a supermassive primordial star forms in a region of the universe with a high 
molecule-dissociating background radiation field, and collapses directly into a 10 4 -10 6 Mq seed BH. 
This also corroborates the results of recent cosmological simulations which suggest that these massive 
BHs were the seeds of a large fraction of the SMBHs residing in the centers of galaxies today. 
Subject headings: Accretion — Black hole physics — Cosmology: theory — early universe — Galaxies: 
quasars: general — Radiation mechanisms: general 



1. INTRODUCTION 

What is the origin of the SMBHs that power the most 
distant quasars? There are currently three main sce- 
narios that are widely discussed for the formation of 
the initial "seed" BHs from which these SMBHs grew 
(e.g. Haiman 2009, 2012; Natarajan 2011; Volonteri 
2010, 2012): (1) massive Population (Pop) III seeds, 
which form from the collapse of ~ 30-300 Mq primor- 
dial stars in dark matter (DM) minihalocs with total 
masses of ~ 10 5 -10 6 M Q at redshifts z > 20 (e.g. Madau 
& Rees 2001); (2) supermassive stellar remnant seeds, 
which form with initial masses of 1O 4 -1O 6 M from the 
direct collapse of ~10 4 K primordial gas in atomic cool- 
ing DM halos with total masses of ~10 7 -10 8 Mq at z 
>10 (e.g. Bromm & Loeb 2003); (3) and seeds formed 
from the collapse of ~T0 3 Mq stars created in runaway 
collisions in dense stellar clusters at z ~ 10-20 (e.g. De- 
vecchi & Volonteri 2009). 

Distinguishing between these three scenarios poses sig- 
nificant challenges, in large part because the high-z 
regime in which these seed BHs are born is too dis- 
tant to be probed directly by existing facilities. How- 
ever, significant progress has been made in detecting 
> 10 9 Mq SMBHs powering quasars at z > 6, the ex- 
istence of which provide significant constraints on the 
nature and growth of their BH seeds (e.g. Willott et al. 
2003; Fan 2006). The strongest constraints to date come 
from the ~ 2 x 10 9 Mq SMBH inferred to be powering a 
quasar at z = 7.085 (Mortlock et al. 2011). Given the < 
800 Myr available for the growth of such a massive BH, 
only the most optimistic models can explain their ori- 
gin from ~100 Mq Pop III seeds, suggesting that more 
massive seeds may have a role to play (e.g. Tyler et al. 
2003; Shapiro 2005; Volonteri & Rees 2006; Tanaka et al. 
2012). 



At the same time, a growing body of theoretical work 
on the supermassive stellar remnant model suggests that 
the conditions required for the formation of supermassive 
stars (SMSs) may be satisfied more often in the early uni- 
verse than previously assumed (Wise et al. 2008; Regan 
& Haehnelt 2009; Shang et al. 2010; Bellovary et al. 
2011; Wolcott-Green et al. 2011; Agarwal et al. 2012; 
Inayoshi & Omukai 2012; Johnson et al. 2012a; Latif et 
al. 2012; Petri et al. 2012). In addition, modeling the 
evolution of rapidly accreting SMSs (e.g. Begelman 2010; 
Hosokawa et al. 2012) and the radiative feedback they 
exert during their growth (Johnson et al. 2012b; but see 
also Dotan & Shaviv 2012) shows that they can attain 
masses > 10 5 Mq before collapsing to BHs. These mas- 
sive BH seeds could in principle grow to SMBHs much 
more quickly than 10 - 10 3 Mq seeds. This is especially 
true in light of numerous other recent theoretical (e.g. 
Volonteri et al. 2005; Pelupessy et al. 2007; Alvarez et 
al. 2009; Milosavljevic et al. 2009a; Noble et al. 2009; 
Park & Ricotti 2011, 2012; Jeon et al. 2012) and obser- 
vational (e.g. Elvis et al. 2002; Wang et al. 2006, 2009; 
Davis & Laor 2011; Bambi et al. 2012; Li et al. 2012) re- 
sults which suggest that accretion onto BHs in the early 
universe was suppressed due to radiative feedback. 

Here we examine the limits that can be placed on the 
initial seed masses and growth of the most massive high- 
z BHs. In the next Section, we survey the possible range 
of time-averaged growth rates and radiative efficiencies 
of high-z SMBHs, in the context of the main models of 
seed formation. In Section 3, we discuss how radiative 
feedback from both accreting BHs and their progenitor 
stars acts to limit the rate of growth of high-z BHs and in 
particular hinders the growth of the least massive BHs. 
Finally, we discuss the implications for the main models 
of seed formation in Section 4. 
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Fig. 1. — The time-averaged Eddington fraction (/Edd/duty) and 
radiative efficiency (e) required for various BH seeds to produce a 
space density ngMBH of SMBHs with masses MBH,final> 10 9 Mq 
by redshift z. Contours show the minimum values of /Edd/duty (1 — 
e)/e required, for three initial BH seeds: M B H,mit=10 5 Mq SMS 
progenitors in M ha i o ,mit=10 7 Mq DM halos (top), 10 2 Mq Pop III 
progenitors in 10 5 Mq halos (middle), and the same case but with 
accretion delayed by 100 Myr due to radiative feedback (bottom). 
The yellow, green and red contours correspond to e=0.07, 0.1 and 
0.15, respectively, for /Edd/duty = 1; super-Eddington growth is 
required in the gray-shaded regions, for these three cases. The 
black crosses denote values of ngMBH f° r SMBHs with masses > 
10 9 Mq inferred from observations. A seed BH with A/BH,init=10 5 
Mq (top) has to accrete continuously (/duty=l) at the Eddington 
rate (/Edd=l) with a radiative efficiency e < 0.14 to explain the 
existence of at least the one z ~ 7 quasar. 



2. OBSERVATIONAL CONSTRAINTS ON 
SUPERMASSIVE BLACK HOLE SEEDS 

We begin by highlighting the observational constraints 
on the accretion history of high-z SMBHs that can be 
derived from the evolution of the comoving number den- 
sity of SMBHs over cosmic time. In particular, we place 
constraints on the time-averaged accretion rates and on 
the efficiency with which radiation is produced during 
accretion, as described below. 

In general, the radiative luminosity, L, of a black hole 
accreting at a rate Mbh is given by L = eMbrc 2 / (1 — e), 
where c is the speed of light and e is the radiative effi- 
ciency, defined as the fraction of the rest mass energy of 
infalling matter that is converted into radiation during 
accretion[3 We express the radiative luminosity as a frac- 
tion /Edd of the Eddington luminosity of the BH, L = 
/Edd^Edd, where i E dd = 1.2 x I0 38 ergs" 1 (M BH /M Q ). 
The accretion rate at which a BH will radiate at a given 
fraction /Edd of its Eddington luminosity ^Edd is then 
given by 

(1 — e)/Edd^Edd 



M, 



BH 



We integrate equation (1) over time to find the final 
SMBH mass -Mbh, final ; &s a function of its initial seed 
mass, MsHjinit; the duty cycle /duty at which it accretes 
(defined as the fraction of time spent accreting), and the 
radiative efficiency e. This yields 



M B H,final = M B H,init X 

/Edd/duty(l ~ e) / ^final ~ fr] 
£ V tedd 



(2) 



exp 



Here £sdd = 450 MyiQ, while i nnal and t init are the ages 
of the universe when the BH attains its final mass and 
at the time of seed formation, respectively. 

With /Edd = 1 and /duty = 1, equation (2) reduces 
to the standard equation describing BH growth at the 
Eddington rate (e.g., Volonteri & Rees 2006). However, 
in Fig. I we have left these two parameters, as well as 
the radiative efficiency, e, free. Instead, we have plotted 
contours of the product /Edd/duty (1 — This allows 

us to constrain this product from the observationally- 
derived number density of BHs. The contours in Fig.l 
show the minimum values that this product must take 
in order for a number density ngMBH of SMBHs with 
masses Mbh, final > 10 9 M Q to be produced by redshift 
z. 

For our comparison with the observational data in 
Fig. 1, we have expressed t nn ai as a function of redshift, 
using the WMAP7 results (Komatsu et al. 2011) for 
the relevant cosmological parameters (J1m = 0.27, 17a = 
0.73 and Hq = 70.3 km s^ 1 Mpc" 1 ), and following the 

1 We note that only a fraction (1-e) of the rest mass of the accret- 
ing material is finally accreted by the BH, as seen by an observer 
at infinity (e.g. Salpctcr 1964; Thorne 1974). The radiation gener- 
ated during accretion which escapes to infinity accounts for the rest 
mass that such an observer concludes is lost due to gravitational 
rcdshifting and time dilation as material falls into the potential 
well of the BH. 

2 This timescalc follows directly from the definition of the Ed- 
dington luminosity: *Edd = CT T C / '47rGm p , where (Tt is the Thom- 
son cross section for electron scattering, G is Newton's constant, c 
is the speed of light and m p is the mass of the proton. 
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approach in Barkana & Loeb (2001). In order to deter- 
mine tinit in equation (2), as a function of tismbh, wc 
have used the Warren et al. (2006) DM halo mass func- 
tion to find the highest redshift at which the comoving 
number density of the DM halos with masses > -Mh a io,init 
which can host BH seed formation is > tismbh- This im- 
plies a lower limit to ii n it , since in principle some of these 
halos could have merged with one another before ifi na i, 
in which case the initial seeds of some of the SMBHs 
would have formed at later times. While we ignore the 
effect of mergers in our modeling, in Section 3.1 we argue 
that mergers play a minor role in the growth of high-z 
SMBHs. Finally, we note that the Warren et al. (2006) 
mass function we have adopted provides a reasonable fit 
at the high redshifts (e.g. z ~ 30) of seed formation to the 
halo mass functions found from cosmological simulations 
by Reed et al. (2007). That said, Fig. 1 also shows that 
our results are relatively insensitive to the exact value 
of fisMBHi except at the highest redshifts where there is 
the least time for SMBHs to form in large numbers (as 
shown by the downturn of the contours at the highest 
number densities). 

In Fig. 1 wc compare these theoretical curves with the 
observational constraints on the evolution of the comov- 
ing number density tismbh of SMBHs over cosmic time. 
We show lower limits on jt-smbh for BHs with masses > 
10 9 Mq at the highest redshifts, as found from detec- 
tions of z ~ 6 quasars in the Sloan Digital Sky Survey 
(SDSS; e.g. Fan et al. 2006, Tanaka & Haiman 2009; 
see also Willott et al. 2010) and from the single quasar 
at z ~ 7 detected in the United Kingdom Infrared Deep 
Sky Survey (UKIDSS; Mortlock et al. 2011)0 At lower 
redshifts, Fig. 1 also shows the number densities of > 10 9 
Mq BHs obtained by integrating the BH mass function 
inferred from the observed luminosity function of active 
galactic nuclei by Shankar et al. (2009). 

In each of the three panels of Fig. 1 the observed num- 
ber densities of SMBHs are plotted together with the 
theoretical contours of /Edd/duty (Te)/e, for three differ- 
ent models of the initial BH seeds and their early evolu- 
tion0 In the top panel, we take values of the initial seed 
mass MeH.init = 10 5 Mq and initial DM host halo mass 
-Mhaio.init = 10 7 Mq, corresponding to the SMS remnant 
model. In the middle panel, we take MsH.init = 10 2 Mq 
and Mhaio.init = 10 5 Mq, corresponding to the Pop III 
remnant model. Finally, in the bottom panel we take the 
same values as in the middle panel but we assume accre- 
tion onto the seed BH is delayed for the first 100 Myr; as 
discussed in Section 3.3, numerous cosmological simula- 
tions suggest that this may be the case for BHs born in 
minihalos at high-z due to radiative feedback from the 
progenitor star and from the accretion process (e.g. Ki- 
tayama et al. 2004; Whalen et al. 2004; Yoshida 2006; 
Alvarez et al. 2009). 

For a given value of e, the contours show the range 

3 We have estimated a conservative lower limit for riSMBH at z ~ 
7 by noting that just one > 10 9 Mq BH was found in the UKIDSS 
Large Area Survey which covered ~ one tenth of the sky, and by 
assuming that such a SMBH could have been detected out to z = 
10. While this is a rough estimate, Fig. 1 shows that our results 
arc relatively insensitive to values of 'ftSMBH in this range. 

4 For simplicity, we consider the two models described in Sec- 
tion 1 which bracket the range of expected initial seed masses: the 
Pop III remnant model and the SMS remnant model. 



of redshifts over which there is sufficient time to grow a 
BH from its initial seed mass to A/bh, final > 10 9 Mq 
assuming accretion at a time-averaged fraction of the 
Eddington rate /Edd/duty In particular, for constant 
Eddington- limited accretion, /Edd/duty = 1, the yellow, 
green and red contours in Fig. 1 correspond to efficien- 
cies e = 0.07, 0.1 and 0.15, respectively. The gray-shaded 
regions show the range of redshifts at which > 10 9 Mq 
BHs must have grown at time-averaged rates exceeding 
the Eddington limit for each of these efficiencies. As 
these gray-shaded regions extend to lower redshifts in 
the Pop III remnant cases (bottom two panels), it is 
clear that relatively low-mass Pop III BHs could only 
be the seeds of the SMBHs inferred at z > 6 in the SDSS 
and UKIDSS if they grew at supcr-Eddington rates for 
a significant fraction of time and/or if their radiative ef- 
ficiency of accretion is relatively low, e < 0.07. However, 
the much more massive SMS remnant seeds could ac- 
crete at sub-Eddington rates and still grow to Mbh, final 
> 10 9 Mq sufficiently rapidly, even allowing for higher 
radiative efficiencies. 

In the next two Sections we consider the values ex- 
pected for the parameters e, /Edd an d /duty, an d what 
they imply for models of SMBH seed formation. 

3. SUPPRESSION OF BLACK HOLE GROWTH VIA 
RADIATIVE FEEDBACK 

Here we discuss various ways in which radiation limits 
the rate of growth of BHs in the early universe, and what 
values are expected for the radiative efficiency e and the 
time-averaged accretion rates /Edd/duty of early BHs. We 
then consider what these limits, taken together, imply for 
the initial seed masses of SMBHs. 

3.1. Limited Growth via Mergers and Lyman- Werner 

Feedback 

In the previous Section we described our model of BH 
growth from a single seed, which neglected growth due to 
mergers. Here we explain why the role of mergers must 
be limited, in part due to radiative feedback. 

Many previous works have examined the possibility of 
SMBHs growing from Pop III seeds via mergers and gas 
accretion (Menou et al. 2001; Haehnelt 2003; Islam et 
al. 2003; Yoo & Miralda-Escude 2004; Volonteri et al. 
2003, 2005; Tanaka & Haiman 2009). If sufficient mass 
can be locked up in the seeds, merging the BHs is an ef- 
fective means of growing a single SMBH that avoids the 
radiative feedback limitations associated with accreting 
gas. Most recently, Tanaka et al. (2012) have reported 
that mergers of Pop III remnant seeds born in DM halos 
with virial temperatures T v i r = 400 K enhance the final 
masses of SMBHs by a factor of 10-1000 However, ac- 
counting for the substantially more massive halos, with 
T vir ~ 1000-2000 K0 in which Pop III star formation is 
found to occur in cosmological simulations (e.g. Yoshida 
et al. 2003; O'Shea & Norman 2007), even at the earliest 

5 These authors also account for the reduced rate of gas accretion 
onto BHs when they are kicked out into the low density regions of 
their host halos due to the emission of gravitational waves during 
mergers (e.g. Favata et al. 2004; Merritt et al. 2004; see also 
Madau et al. 2004). 

6 The virial temperature of a halo increases with halo mass as 

2 

T vir <x M h 3 alo (e.g. Barkana & Loeb 2001). 
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epochs (e.g. Gao ct al. 2007), implies a decrease in the 
host halo abundance by a factor of 1-2 orders of mag- 
nitude (as found from large-scale cosmological simula- 
tions of hierarchical structure formation; e.g. Reed et al. 
2007). Assuming the same star formation efficiency and 
Pop III initial mass function from Tanaka ct al. (2012) 
this, in turn, suggests that the number of mergers that 
occur in the assembly of a high-z SMBH is 1-2 orders of 
magnitude lower than Tanaka et al. find, which implies 
that mergers are likely to contribute only minimally to 
the growth of SMBHs at high-z (see also Madau et al. 
2004) . We also note that the results of recent large-scale 
cosmological simulations tracking the build-up of SMBHs 
at high redshift also suggest that only a negligible frac- 
tion of their mass is acquired in mergers (DcGraf et al. 
2012). 

The rate of mergers that grow SMBHs from Pop III 
seeds is further reduced due to the build-up of the 
molecule-dissociating Lyman-Werner (LW) radiation 
field which acts to slow the collapse of primordial gas 
in DM halos, and thus lowers the Pop III star formation 
rate (SFR) (e.g. Haiman et al. 1997; Glover & Brand 
2001; Machacek et al. 2001; Ricotti ct al. 2001; Ciardi 
& Ferrara 2005; Mesinger et al. 2006; Wise & Abel 2007; 
O'Shca & Norman 2008). The most recent cosmological 
simulations tracking the build-up of the LW radiation 
field suggest that the Pop III SFR is reduced by a fac- 
tor of a few compared to the SFR in the absence of LW 
feedback (Aim et al. 2012; Hummel et al. 2012; Johnson 
et al. 2012a; Wise et al. 2012). Accounting for this even 
further reduction in the number density of Pop III seed 
BHs, we can safely conclude that mergers are likely to 
be responsible for only a small portion of the growth of 
SMBHs. The vast majority of the growth must be due 
instead to accretion of gas (and perhaps to a much lesser 
extent DM; e.g. Hu et al. 2006, Guzman & Lora-Clavijo 
2011). This justifies our simplified approach of focus- 
ing on the growth of BHs due to accretion only, and not 
incorporating the effects of mergers. 

A lower rate of (gas poor) mergers of seed BHs also 
suggests higher values for the spins of the seeds (e.g. 
Gammic et al. 2004; Volonteri et al. 2005, 2012; Berti 
<fc Volonteri 2008), since mergers of BHs with randomly 
oriented spins will tend to spin down fast rotating BHs. 
This translates into a higher radiative efficiency of accre- 
tion, since more energy can be extracted from the spin 
of the BH (Blandford & Znajek 1977) and from the hot- 
ter accretion disk that extends further inward towards 
the horizon of a faster-spinning BH (e.g. Novikov & 
Thorne 1973) Q Next, we discuss the expected radiative 
efficiency of high-z SMBHs and the impact this has on 
their growth. 

3.2. Low Eddington Accretion Rate due to High 
Radiative Efficiency 

Following equation (1), for higher e the Eddington lu- 
minosity is generated at lower accretion rates, which 
implies that Eddington-limited accretion proceeds at a 

7 It is worth noting that the spins of local SMBHs arc found 
to be > 60 percent of the maximum allowed value (Brenneman et 
al. 2011), which suggests that their radiative efficiencies are higher 
than the value expected for a non-rotating BH, e > 0.06 (e.g. Noble 
et al. 2011). 



lower rate for higher radiative efficiencies. For higher e 
a BH must accrete at a higher time-averaged accretion 
rate /Edd/duty to grow to a given mass within a given 
time, as shown in Fig. 1. 

Given that the Eddington accretion (defined as that 
which produces the Eddington luminosity) rate is sen- 
sitively dependent on e and that the radiative efficien- 
cies of BHs can range from ~ 0.025 to ~ 0.4 (see e.g. 
Milosavljevic et al. 2009b and references therein)0, it is 
vital to constrain this quantity in order to understand 
the growth of the earliest SMBHs (e.g. Shapiro 2005; 
King & Pringle 2006; Volonteri & Rees 2006; Tanaka et 
al. 2012). Fortunately there is a growing body of ob- 
servational evidence which provides some guidance, and 
in particular suggests that SMBHs at high-z tend to be 
rapidly rotating and have relatively high radiative effi- 
ciencies. Numerous authors have argued that in order 
to account for the observations, SMBHs must have typ- 
ical values of e > 0.1-0.15 (e.g. Elvis ct al. 2002; Yu 
& Trcmaine 2002; Volonteri ct al. 2005; Shankar et al. 
2010; but see e.g. Raimundo et al. 2012), and that e 
tends to increase with redshift (Wang et al. 2006, 2009; 
Barausse 2012; Volonteri et al. 2012; see also Maio et 
al. 2012) and with BH mass (e.g. Davis & Laor 2011; 
Shankar et al. 2011), with inferred values as high as 
e - 0.3-0.4 for > 10 9 M Q BHs at high-z (Li et al. 2012). 

As shown in Fig. 1, for such high efficiencies (e.g. e 
> 0.15) ^100 Mq Pop III seed remnants would have to 
grow at super-Eddington time-averaged accretion rates 
(i.e. with /Edd/duty > 1) in order to explain the Mort- 
lock et al. (2011) SMBH, the SMBHs powering the SDSS 
quasars and other >10 9 Mq BHs at z > 4 (to have 
/Edd/duty(l — e ) / e larger than the lower limit shown in 
Fig. 1). However, more massive (> 10 5 Mq) SMS rem- 
nants could grow sufficiently fast to explain these high-z 
SMBHs, even at sub-Eddington rates. 

Fig. 1 also shows that for the highest radiative effi- 
ciencies e ~ 0.3-0.4, even the supermassive seeds from 
SMS would have to accrete at super-Eddington time- 
averaged rates (i.e. /Edd/duty > 1) just to grow to 
M BH ,finai > 1O 9 M by z ~ 2-4. This suggests that 
SMBHs with such high radiative efficiencies at higher 
rcdshifts may have undergone a period of rapid accre- 
tion with a lower radiative efficiency in the past. While 
we cannot know the entire accretion history of a given 
SMBH at high-z, we can place some constraints on 
the amount of high energy radiation emitted during its 
growth from the observationally-infcrred size of the H II 
region surrounding it. For the case of the z ~ 7 quasar, 
Mortlock et al. (2011) report that it resides in an H II 
region that is ~ 1.9 physical Mpc in extent, which is 
smaller than those typically found surrounding quasars 
at z > 6. This suggests that accretion may indeed have 
been radiatively inefficient during much of the growth of 
the BH. Alternatively, however, the accreting BH could 
have instead been obscured for much of its lifetime, re- 
sulting in a small fraction of ionizing photons escaping 
into the intergalactic medium (see Bolton et al. 2011). 
This would be consistent with the merger-driven model 

8 These values are expected for accretion through geometrically 
thin disks. For spherically symmetric or advection dominated ac- 
cretion (as expected at low accretion rates; e.g. Narayan & Yi 
1995), the radiative efficiency can be much lower. 
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for SMBH growth (e.g. Sanders et al. 1988; Hopkins ct 
al. 2006) as well as with recent observational evidence 
suggesting that a large fraction of accreting SMBHs at 
high-z are in fact buried within significant amounts of 
gas and dust that prevent the escape of ionizing radia- 
tion (e.g. Fiore et al. 2012; see also Kelly & Shen 2012). 

3.3. Sub-Eddington Accretion due to Progenitor- and 
Accretion- Generated Radiation 

By definition, accretion may in principle be limited 
to the Eddington rate due to electron scattering of the 
photons generated in the accretion process, but other 
radiative processes can further limit the accretion rate. 
Milosavljevic et al. (2009a, b) find that photoheating and 
pressure from ionizing and Lyman-a photons render the 
accretion of gas onto a BH intermittent (/duty < 1), with 
a time-averaged accretion rate of just ~ 0.3 times the 
Eddington rate (i.e. /Edd/duty ~ 0.3). Park & Ricotti 
(2011, 2012) report similar results, although they find 
that accretion at the Eddington rate can be achieved for 
sufficiently dense accreting gas (see also Li 2011). The 
results of larger-scale cosmological simulations also sup- 
port the conclusion that accretion onto early BHs occurs 
at sub-Eddington rates due to accretion-generated radia- 
tive feedback (e.g. Pelupessy et al. 2007; Alvarez et al. 
2009; DeGraf et al. 2012; Jeon et al. 2012). 

An additional bottleneck to efficient accretion onto 
Pop III seeds born in minihalos is that the intense ion- 
izing radiation emitted by their progenitor stars drives 
dense gas out of the halo (Kitayama et al. 2004; Whalen 
et al. 2004; Alvarez et al. 2006), leaving the BH in a 
low-density medium from which it cannot accrete rapidly 
(Yoshida 2006; Abel et al. 2007; Johnson & Bromm 
2007)0 This results in accretion rates orders of mag- 
nitude below the Eddington rate for up to ~ 10 8 years 
before dense gas recollapses into the halo. The bottom 
panel of Fig. 1 shows the effect of such a delay in accre- 
tion onto a 100 Mq Pop III seed BH initially formed in 
a 10 5 Mq DM halo. Comparing this to the case without 
a delay (middle panel) shows that the effect is compara- 
ble to a decrease of ~ 10-20% in the average accretion 
rate onto the seeds of SMBHs formed by z ~ 6-8. Thus, 
the radiative feedback from Pop III progenitor stars can 
significantly slow down the growth of Pop III seeds to 
SMBHs. 

While radiative feedback from both accretion and the 
progenitor star are likely to limit Mbh to sub-Eddington 
values at early times, we do note that SMBH-powcred 
quasars at z ~ 6 are inferred to have Eddington ratios 
(/Edd) near unity, although their duty cycle (/duty) is not 
well-constrained (Willott ct al. 2010; see also Shankar et 
al. 2010). Overall, however, we conclude that radiative 
feedback, at least at the earliest times, appears likely to 
keep the time-averaged accretion rate of Pop III remnant 
seeds at / E dd/dut y < 1. 

Accretion-generated radiative feedback may also limit 
the growth of SMS remnant BH seeds to sub-Eddington 
rates (e.g. Johnson et al. 2011). However, as the in- 
tense ionizing radiation emitted by rapidly growing SMSs 

9 Less massive (20—40 Mq) Pop III progenitor stars may not 
evacuate the gas as completely, but the BHs they create are also 
likely to be ejected from their host halos due to kicks they receive 
during core collapse (Whalen & Fryer 2012). 
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Fig. 2. — The minimum initial BH seed mass MBH,init required to 
form a JV/bh, final > 10 9 Mq SMBH by rcdshift z, at time -averaged 
Eddington fractions of /Edd/duty = 1 {solid lines) and 0.5 (dashed 
lines), for three different radiative efficiencies: e = 0.15 (red), 0.1 
(green), and 0.07 (yellow). The three panels correspond to the 
cases shown in the panels of Fig. 1. The largest seed masses 
are implied for the case of seed formation in a A'/ na i ; n j t =10 5 
Mq halo and accretion delayed by 100 Myr due to radiative feed- 
back (bottom panel), while less stringent constraints are given for 
the cases with no delay (middle panel) and with larger initial 
halo mass (M nal0jinit = 10 7 Mq instead of M halo , init = 1O 5 M ) 
(top panel). For the high radiative efficiencies inferred for high-z 
SMBHs, these objects likely grew via constant (or time-averaged) 
supcr-Eddington accretion and/or started as massive Mbh init^ 
10 3 -10 5 Mq BH seeds. 
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cannot escape their heavy accretion flows, their host ha- 
los are less likely to be photoevaporated (Johnson et al. 
2012b; see also Hosokawa et al. 2012); in turn, the BH 
remnants they leave behind are likely to accrete more 
rapidly than Pop III remnant seeds. This is in basic 
agreement with the suggestion by Salvaterra et al. (2012) 
that more massive seed BHs accrete with higher Edding- 
ton fractions than do lower mass seeds. 

4. IMPLICATIONS AND DISCUSSION 

In the previous Sections we have highlighted the con- 
straints that observations of high-z quasars place on the 
nature of the BH seeds and subsequent growth of SMBHs 
in the early universe, and we have reviewed the ways in 
which radiative feedback from stars and accreting BHs 
is expected to limit the growth of these objects. We now 
discuss the implications of these findings for the main 
models of SMBH seed formation. 

In Section 3.1 we argued that mergers of BH seeds 
are likely to play only a minor role in growing SMBHs, 
and that the rate of accretion required to grow the ob- 
served SMBHs (i.e. the quantity /Edd/duty (l-e)/e) de- 
pends weakly on their number density (see Fig. 1). With 
these assumptions we take a simplified approach and 
solve equation (2) for the initial seed mass A/BH.init as 
a function of e, /Edd, and /duty Figure 2 shows the min- 
imum BH seed mass, MBH,init, required to grow a SMBH 
with a mass > 10 9 Mq by redshift z, with the three pan- 
els corresponding to the scenarios shown in the panels 
in Fig. 1. In each panel, the minimum BH seed mass 
is shown for two choices of the time-averaged Eddington 
fraction (/Edd/duty = 0.5 and 1) and for the same three 
radiative efficiencies highlighted in Fig. 1 (e = 0.07, 0.1, 
and 0.15). For each of these cases we have assumed a 
space density of SMBHs of tismbh = 1 Gpc -3 (comov- 
ing), but as Fig. 1 shows the results are not strongly 
sensitive to this choice. 

The top two panels differ by very little, which is a 
reflection of the fact that the time between the for- 
mation of the first 10 5 Mq halo and that of the first 
10 7 M Q halo in a 1 Gpc 3 comoving cosmological vol- 
ume is small compared to the time from their forma- 
tion to redshift z (on the horizontal axis). The values 
of MBH.init shown in the bottom panel, however, are 
somewhat higher, reflecting the fact that the 100 Myr 
delay assumed in this case is longer than the time be- 
tween the formation of the first 10 5 and 10 7 Mq halos. 
In every case, it is clear that high radiative efficiencies 
(e > 0.1-0.15) such as those inferred from observations 
of high-z SMBHs (see Section 3.2) imply seed masses for 
observed SMBHs at z > 7 of M B H,init > 10 3 -10 5 M , as- 
suming constant Eddington-limited accretion (/Edd/duty 
= 1). For lower time-averaged accretion rates, as sug- 
gested by simulations of BH accretion (see Section 3.3), 
the implied minimum seed masses are much higher, up 
to M B H,init > 10 6 M Q for /Edd/duty ~ 0.5 and over the 
same range of radiative efficiencies. 

If accretion is radiatively efficient and the time- 
averaged accretion rate is sub-Eddington, as discussed in 
Section 3, then the initial seed masses of the highest red- 
shift SMBHs must have been very high, perhaps exceed- 
ing even the range predicted by the SMS remnant model 
(10 4 -10 6 M Q ) in order to explain the 2 x 10 9 M Q BH 



reported at z ~ 7 (Mortlock et al. 2011). Interestingly, 
this is also consistent with the results of recent large-scale 
cosmological simulations tracking SMBH growth which 
suggest that such high-z quasars can be explained by 
starting with initial seed BH masses of 10 5 M© (Li et 
al. 2007; Di Matteo et al. 2012). Taken together, we 
conclude that the available theoretical and observational 
evidence strongly suggests that the seeds of SMBHs at 
high-z were, at least in some cases, very massive. In turn, 
as the SMS remnant model produces the most massive 
seeds of the three main models discussed in the introduc- 
tion, the evidence suggests that this model may be the 
most viable (see also e.g. Natarajan & Volonteri 2012). 

One of the arguments against the SMS remnant model 
has been that SMSs are exotic objects, and they may 
never form in our universe. A growing body of work 
is elucidating the conditions required for these objects 
to form (e.g. Bromm & Locb 2003; Koushiappas et al. 
2004; Lodato & Natarajan 2006; Begelman et al. 2006; 
Spaans & Silk 2006; Regan & Haehnelt 2009; Shang et al. 
2010; Schleicher et al. 2010; Ball et al. 2011; Wolcott- 
Green et al. 2011; see also Mayer et al. 2010 on massive 
seed formation from metal-enriched gas), and there are 
suggestions that the conditions for SMS formation may 
occur much more often in the early universe than previ- 
ously assumed (Dijkstra et al. 2008; Agarwal et al. 2012; 
Johnson et al. 2012a; Petri et al. 2012). Indeed, Agar- 
wal et al. (2012) find that a sufficient number of SMSs 
may form to provide the seeds for a large fraction of the 
SMBHs in the centers of galaxies today. This develop- 
ment offers a completely independent reason to seriously 
consider SMS remnants. 

While the SMS remnant model offers an explanation 
for the origins of the earliest SMBHs that is consistent 
with available data, we note that other models start- 
ing with lower initial seed masses can, in principle, also 
explain the data. For instance, as shown in Fig. 2, if 
A^BH.init ~ 10 2 Mq seeds accrete gas with very low radia- 
tive efficiency (e < 0.07) continuously at the Eddington 
rate (or, perhaps intermittently, above it; e.g. Jaroszyn- 
skietal. 1980; Collin et al. 2002; Kawaguchi et al. 2004; 
Ohsuga et al. 2005; Volonteri & Rees 2005; Kurosawa & 
Proga 2009; Wyithe & Loeb 2011; Begelman 2012a; Li 
2012) then, even allowing for a ~ 100 Myr delay due to 
radiative feedback (bottom panel) , they could grow to > 
10 9 Mq by z > 7. However, it is possible that a large 
fraction of the mass in such super-Eddington flows is lost 
to a wind instead of being accreted onto the BH (see e.g. 
Begelman 2012b; Dotan & Shaviv 2011). It also remains 
to be demonstrated that these conditions are likely to be 
realized at z > 7. 

It is difficult to verify the SMS remnant model without 
direct observational evidence of the existence of SMSs. 
We note, however, that there are observational signa- 
tures of SMSs that may be detected by future missions 
such as the James Webb Space Telescope (e.g. Johnson 
et al. 2012b). They may also leave unique chemical sig- 
natures that could be detected in Lyman-limit systems 
(Woosley 1977; Fuller & Shi 1997). In addition, upon 
their collapse SMSs are predicted to emit a large neu- 
trino flux that could be detectable (by e.g. IceCube) (Shi 
& Fuller 1998; Linke et al. 2001; Fryer & Heger 2011; 
Montero et al. 2012), as well as producing gravitational 
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wave signatures that could be uncovered by the Laser In- 
terferometer Gravitational Wave Observatory (e.g. Fryer 
& New 2011)0 and extremely bright supernovae that 
could be found by the Wide-Field Infrared Survey Tele- 
scope (e.g. Whalen ct al. 2012). In lieu of such ob- 
servations, the best evidence for the existence of SMSs 
remains the SMBHs to which their remnants may have 
grown at the highest rcdshifts. 
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